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Abstract 

The effects of shell closure in nuclei via the cluster decay is studied. In 
this context, we have made use of the Preformed Cluster Model (PCM) of 
Gupta and collaborators based on the Quantum Mechanical Fragmentation 
Theory. The key point in the cluster radioactivity is that it involves the 
interplay of close shell effects of parent and daughter. Small half life for 
a parent indicates shell stabilized daughter and long half life indicates the 
stability of the parent against the decay. In the cluster decay of trans lead 
nuclei observed so far, the end product is doubly magic lead or its neighbors. 
With this in our mind we have extended the idea of cluster radioactivity. We 
investigated decay of different nuclei where Zirconium is always taken as a 
daughter nucleus, which is very well known deformed nucleus. The branching 
ratio of cluster decay and a-decay is also studied for various nuclei, leading 
to magic or almost doubly magic daughter nuclei. The calculated cluster 
decay half-life are in well agreement with the observed data. First time a 
possibihty of cluster decay in "^^^U nucleus is predicted. 
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1 Introduction 



When a charged particle heavier than e, but hghter than a fission fragment 
is emitted by an unstable nucleus, the process is called cluster radioactivity 
or heavy ion radioactivity. Cluster radioactivity was predicted theoretically 
in 1980 by Sandulescu, Poenaru and Greiner [1] on the basis of Quantum 
Mechanical Fragmentation Theory (QMFT). The first experimental identi- 
fication of a case of radioactive decay of heavy nuclei by emission of a nuclear 
fragment heavier than was done by Rose and Jones [2], who studied the 
decay of ^^^i?a by ^^C emission with a half-life of (3.7 ±1.1) x 10^ yr. Rose 
and Jones measured a branching ratio of B(^'^C/a) = 8.5(2.3) x 10~^° . This 
result was then confirmed by Alekshandrov et al. [3], Gale et al. [4] and 
Price et al. [5]. 

Cluster radioactivity has grown and there arc more than 24 cases of clus- 
ter radioactivity with partial half-lives ranging from IoqiqT = > 3.63(s) to 
> 29.20(s) [6]. There exists a whole family of cluster decay modes like ^'^C 
radioactivity, ^'We radioactivity, "^^Mg radioactivity and so on. The mother 
nuclei range from '^'^^Fr to ^^^Cm, all from trans-lead region and cluster ra- 
dioactivity in this region indicates the presence of heavier clusters. Shell 
effects arc clearly manifested in cluster radioactivity since experiments show 
that the cluster decay result in magic daughter products. 

The study of a-decay and cluster decay has been used for identifying the 
shell closure effects including even the very weak sub-shell closures [7, 8, 9, 
10]. Both the cases of large and small decay half-lives are important. Small 
ones imply closed shell effects in daughter nucleus and the large ones indi- 
cate closed shell effects of the mother nucleus. In decay calculations presence 
of the known spherical or deformed daughter should result in a small decay 
half-life or a small decay half-life should refer to the existence of a known or 
new, spherical or deformed, closed shell for the daughter nucleus. 

Superdeformation in medium-mass nuclei is a subject of interest and it 
was first discovered in the actinide fission isomers [11]. Later it was explained 
from the secondary minimum at very large deformation[12]. 7- ray transtions 
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between states in the N — Z nucleus ^°Zr was studied to see the deforma- 
tion, and it was observed that this is the most deformed nucleus in nature, 
with a quadrupole deformation ^2 ~ 0.4 [13]. Superdeformation band in the 
N,Z — AO mass region has also been studied experimentally by the Baktash 
et.al., [14] and a significant subshell closure at the = 40 neutron number 
has been studied in the ^^Ni [15]. 

We have not performed the cluster decay calculations, only to search the 
spherical and/or deformed magicity i.e. not to just study the stability /insta- 
bility of the concerned nuclei, but also for looking the most probable cluster 
decay modes of ^^^U nucleus that has recently been studied for a-decay [16]. 

We are using the Preformed Cluster Model(PCM) of Gupta and Collab- 
orators, where the cluster is assumed to be preformed in the mother nucleus. 
In this model the preformation probability (also known as spectroscopic fac- 
tor in various models) for different possible clusters are calculated by solving 
the Schrodinger equation for the dynamical flow of mass and charge.Another 
such model based on the calculations of spectroscopic factors for different 
clusters, has also been given Blendowske and Walliser [17, 18]. In this model 
they have given parameterization of the spectroscopic factor and meanwhile, 
other contributions have extended these results [19]. Also, the so called " 
Semimicroscopic Algebraic Cluster Model" [20, 21, 22, 23] contributes to pa- 
rameterization of the spectroscopic factors. Deformations of the neutron rich 
clusters have effect on penetration probability as the inclusion of the deforma- 
tions reduces the barrier height. This fact is incorporated into PCM through 
parameter R, which results in the lowering of barrier height [24]. Recently 
[25, 26], role of deformation on binary and ternary clusterization has been 
studied to investigate the exotic nuclear shapes, on the basis of U(3)selection 
rule . 

The paper is organised as follows. The calculations are made by using the 
preformed cluster-decay model (PCM) of Gupta and collaborators [24, 27, 28] 
whose outhne is presented in section II. Section III deals with the calculations 
and results obtained from this study. 



3 



2 The Preformed Cluster Model (PCM) 



The Preformed Cluster Decay Model (PCM) [24, 27, 28] is based on the Quan- 
tum Mechanical Fragmentation Theory [29, 30, 31]. The dynamical collective 
coordinates of mass rj and charge asymmetries rjz are in addition to the usual 
coordinates of relative separation R. 

The decay constant A (or the decay half-life T1/2 ) in PCM is defined as 

A = 1^ = Po^oP. (1) 
-i 1/2 

Here Pq is the cluster (and daughter) preformation probability and P the 
barrier penetrability which refer, respectively, to the rj and R motions. The 
z/q is the barrier assault frequency. The Pq are the solutions of the stationary 
Schrodinger equation in 77, 

which on proper normalization are given as 

Po^^\^P^'\v{A^))n2/A), (3) 

with i=l or 2 and z/=0, 1,2,3.... Eq. (2) is solved at a fixed R = Ra = Ct{= 
Ci + C2), the first turning point in the WKB integral for penetrability P 
(Eq. 5). The Ci are Siissmann central radii = Ri — (l/Ri), with the radii 
Ri = 1.28^-/^ - 0.76 + O.SA;^^^fm. 

The fragmentation potential Vr(?7) in (2) is calculated simply as the sum 
of the Coulomb interaction, the nuclear proximity potential [32] and the 
ground state binding energies of two nuclei, 

ViRa,v) = EBiA, Zi) + + Vp, (4) 

with B's taken from the 2003 experimental compilation of Audi et al., [33] 
and from the 1995 calculations of MoUer et al. [34] whenever not available in 
[33]. Thus, full shell effects are contained in our calculations that come from 
the experimental binding energies and/or from the calculations of MoUer et 
al. [34]. The charges Zi and Z2 in (4) are fixed by minimizing the potential 



4 



in r]z coordinate with Z. The Coulomb and proximity potentials in (4) are 
for spherical nuclei. The mass parameters ^^^(77), representing the kinetic 
energy part in (2), are the classical hydrodynamical masses [35]. 

The WKB tunnelling probability is P = PiPb with 

P^ = exp[-l r {2i^[V{R) - V{R^)]y/^dR] (5) 

Pb = exp[~ r{2i^[V{R) - Q]Y'^dR]. (6) 
nJ Ri 

These integrals are solved analytically [28] for it!^, the second turning point, 
defined by V{Rb) — Q- value for the ground-state decay. 
The assault frequency vq in (1) is given simply as 

= {2E^/^Y'^/Ro. (7) 

with E2 = {Ai/A)Q, the kinetic energy of the lighter fragment, for the Q- 
value shared between the two products as inverse of their masses. Here Rq 
is the equivalent spherical radius of the mother nucleus 

3 Results and Discussion 

In this paper we have performed three types of calculations. Firstly, we have 
focused on the search for the magicity (spherical and /or deformed) taking 
the Zirconium nucleus as daughter in the decay of mother nuclei, which are 
taken from mass region A = 88 to A = 186. Secondly, we have investi- 
gated the cluster decay of various nuclei taken from the heavy mass region 
and have also compared data with experimental results. Thirdly, we have 
searched for the possible cluster decay modes of "^^^U nucleus, whose alpha 

decay has been studied recently [16], along with few other isotopes of Ura- 
nium (230,232,234,236 

In the first part of our calculation we studied the decay of different nu- 
clei taken from the intermediate mass region. We considered the Zirconium 
nucleus as daughter product in the decay of all these nuclei. Mother of even 
mass number is chosen and it results in a daughter with different even neu- 
tron number. The decaying nuclei are from ^^Ru, ^^-^^Pd,^^-^^Cd,^^-^^^Cd, 
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94-116^^^ ioo-i26j.g^ i58-i66j.g io4-i34^g^ 1^8-176^^ ^ ^''^-^^''Ba. Selection of 
isotopes of different mother nuclei for the calculations has been made on the 
basis of Q value and the decay of only those was considered for which it is 
positive. 

The histogram in Fig.(l), contains maxima and minima of half-lives of 
above mentioned parent nuclei with neutron number of the daughter. The 
daughter was fixed to be Zirconium nucleus i.e. Zi — 40. We considered 
all possible cluster decay channels of above mentioned nuclei. Half hves for 
different decaying parents giving Zi — 40 and corresponding neutron number 
A^i of the daughter with the complementing cluster were calculated. There 
exists different combinations of parents and clusters giving rise to Zi — 40 
with a particular value of A^i i.e. a particular isotope of Zirconium, as a 
daughter product. Among different values of half lives of all such possible 
decays resulting into the given Zirconium isotope, we noted the maximum 
and minimum values only. Then these two extreme values of the half life 
were plotted for different isotopes of Zirconium. 

In another histogram the Fig. (2), in which even neutron number of 
daughter nuclei vs half-lives of various mother nuclei are shown. The height 
of the bar corresponding to the daughter neutron number represents the sta- 
bihty/ instability of the daughter nucleus. In this calculation the shell closure 
effects appear at A^i = 40 and at A^i = 82. Also at A^i = 70, these appear, 
but weakly. 

Shell effects can be seen in Fig. (3) also, where we have displayed half-lives 
for the following decay modes: ^^C, '^^Ne, versus the neutron number 

of the daughter, N^. The half-lives shows the minima at the magic number 
for all the decay mode considered above. The Q-value is one of the physi- 
cal quantity which plays a very important role in any spontaneous nuclear 
decay with emission of charged particles. The Q-value behavior has been 
noticed for these calculations, as the size of the cluster increases Q- value also 
increases [7] . The calculations belonging to the same atomic number of the 
daughter are joined with a line of a style mentioned in the figure for Z^ — 82 
and 84. For the ^^C radioactivity where daughter is always with same atomic 
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number and neutron numbers vary, we have observed in our calculations that 
Nd — 126 appears as a magic number in all the cases, same is for the 
cluster when the daughters are Z^, — 80, 82. A strong shell effect can be seen 
in this figure. As a rule, the shortest value of the half-hfe is obtained when 
the daughter nucleus has a magic number of neutrons N^ — 126 and protons 
Zd = 82. 

Same calculation has been made but for different clusters taking the same 
atomic number of different daughter nuclei Z^ — 80, 81, 82, 83, 84 . For the 
^^O radioactivity where daughter is always with same atomic number and 
neutron number vary, we have observed that Na — 126 appears as a magic 
number in all the cases, same is for the '^'^Ne cluster when the daughter are 
Zd — 80, 81, 82. Again a strong shell effect can be seen in Fig.?,. The most 
probable cluster from heavy nuclei in PCM model and their characteristics 
(i.e., Q-value and Half-life) are compared with the experimental data (given 
in [6, 36, 37, 38] and the references therein) in Table 1. 

In the third part of the calculations, we have predicted for the first time 
the possible cluster decay modes of ^^^C/ which was recently investigated for 
alpha-decay [16]. The a-decay calculation of ^^^C/ in the Quantum Mechani- 
cal Fragmentation Theory has been studied and we have calculated the decay 
properties not only for ct-decay but also for other possible clusters. The half 
life for for a-decay, using PCM calculations comes out to be Ti/2 = 4.36ms 
and its experimental value is Ti/2 — 0.51ms. The calculated Q-value is, 
^PCM ^ 8.788Mey and its experimental value is E^^p = 8.612Mey. Shell 
closure effects play an important role in the cluster decay studies. Predicted 
half-hves 7i/2(s) and other characteristics for this isotope of Uranium are 
given in table(2). ^Be, ^^C and ^^~^^A'"e appear as the most probable cluster 
decay modes in our calculations. 
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Figure Captions 



Fig. 1: Predicted half-lives minimum and maximum with daughter (Zr) 
neutron number for various cluster decay mode. 

Fig. 2: Predicted half-lives (minimum only) plotted as a function of the 
daughter neutron number for various cluster decay mode. 

Fig. 3: Predicted half-lives for the following decay modes: ^^C, ^^"^'^O, 
^^Ne, versus the neutron number of the daughter, Nfi. 
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Table 1: Comparision of experimental Q- values, Half-lives with calculated 
Q-values, Half-lives in PCM a-decay of cluster emitters. 



Parent nucleus 
Z and A 


Emitted cluster 
Zf. and Af. 


Qexp. 

MeV 


QFUM. 

MeV 


rpexp 
-^1/2 

sec. 


rpFCM 

^1/2 

sec 




87 - 221 


6- 14 


31.294 


31.292 


14.53 


18.044 




88 - 221 


6 - 14 


32.402 


32.395 


13.39 


17.25 




88- 222 


6- 14 


33.052 


33.05 


11.01 


15.098 




88- 223 


6-14 


31.839 


31.829 


15.15 


18.548 




88- 224 


6- 14 


30.541 


30.536 


15.69 


19.676 




88 -226 


6- 14 


28.198 


28.197 


21.22 


25.805 




89- 225 


6 -14 


30.479 


30.477 


17.16 


22.664 




90- 226 


8-18 


45.731 


45.726 


>16.76 


23.289 




90- 228 


8- 20 


44.730 


44.724 


20.72 


23.357 




90- 230 


10- 24 


57.765 


57.759 


24.61 


26.849 




90 - 232 


10 - 24 


54.491 




>29.20 






90 - 232 


10 - 26 


55.973 


55.964 


>29.20 


29.831 




91- 231 


10 -24 


60.413 


60.409 


22.88 


24.999 




91 -231 


9-23 


51.854 


51.844 


>26.02 


26.463 




92- 230 


10- 22 


61.390 




19.57 






92- 232 


10- 24 


62.312 


62.309 


20.42 


21.528 




92 -233 


10- 25 


60.736 




24.84 






92- 233 


10- 24 


60.490 




24.84 






92 -233 


12- 28 


74.235 




27.59 






92- 234 


12- 28 


74.118 


74.109 


25.74 


27.746 




92- 234 


10- 24 


58.831 


58.825 


>25.92 


28.775 




92 -234 


10- 26 


59.473 


59.465 


25.92 


26.411 




92- 235 


10 -24 


57.358 


57.362 


27.42 


32.511 




92- 235 


10-25 


57.717 


57.756 


27.42 


31.173 




92- 235 


12- 28 


72.162 




>28.09 






92 - 236 


12- 28 


70.558 




27.58 






92- 236 


12- 30 


72.280 




27.58 






93-237 


12-30 


74.791 


74.816 


>27:57 


28.524 




94- 236 


12-28 


79.674 


79.668 


21.67 


21.637 




94- 238 


14- 32 


91.198 


91.189 


25.27 


27.069 




94- 238 


12- 28 


75.919 


75.91 


25.70 


28.395 




94- 238 


12- 30 


76.801 


76.822 


25.70 


25.542 




94- 240 


14- 34 


91.038 


91.027 


>25:52 


25.904 




95- 241 


14- 34 


93.931 


93.925 


>25:26 


23.204 




96- 242 


14- 34 


96.519 


96.509 


23.15 


21.08 
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Table 2: Predicted half-lives 7i/2(s) and other characteristics for cluster de- 
cays of Uranium nuclei. For the Q- value estimates the masses for these par- 
ents are taken from the Audi et. al. mass table [33] and MoUer- Nix [34]. 



Pa rpnt 


Flmitted 


Dauphter 


J-'<-'gi -'-1/2 


Pfpf nrTTi f ,1 on 

J. ± VyJ-WX XJ-XCu UJ.V./XJ. 


F)pcav 


O V/^lnp 




nucleus 




nucleus 


sec 


Drobabilitv 

-L k_/ <_Ai Kj ILL \j y 

Po 


constant 


MeV 




218U 


"^He 


2i4Th 


-2.36 


2 50x10""^ 


1 59x10^2 


8.788 






^Be 


2iORa 


19.759 


9 72x 10~^6 


1.21x10-2° 


16.521 






i°Be 


Ra 


74.199 


4 68x 10"^^ 


4 38xl0-'''5 


7.604 






12Q 


Rn 


24.992 


1.07x 10"22 


7 06x10-26 


31.037 






ISQ 


200pQ 


41.371 


7 68x 10"^° 


2 95x10-42 


39 658 






20O 


198po 


58.519 


4.09x10-33 


2.10x10-5*^ 


33.6 






22Ne 


196pb 


39.92 


2.99 X 10-32 


8.33x10-41 


55.307 






24Ne 


194pb 


46.19 


2.83 X 10-33 


4.48x10-47 


52.081 






28Mg 


190 Hg 


42.169 


3.85 X 10-34 


4.70x10-43 


68.311 




230 u 




226Th 


8.024 


1.15X 10-°^ 


6.56X 10-°9 


5.993 






lOBe 


220Ra 


64.1 


8.55X 10^22 


5.50x10-65 


8.736 






14(3 


2i6Rn 


30.109 


3.04X 10-22 


5.40x10-31 


28.34 






20O 


210Th 


28.69 


2.06X 10-23 


1.42x10-29 


43.771 






23 p 


207Bi 


34.014 


1.48X 10-25 


6.71x10-35 


48.339 






24Ne 


206p]^ 


24.286 


7.45 X 10-24 


3.59x10-25 


61.35 








204p^ 


33.291 


1.62X 10-26 


3.55x10-34 


56.294 




232 y 


^He 


228Th 


10.913 


5.41 X 10-°9 


8.47x 10-12 


5.414 






lOBe 


222Ra 


74.01 


4.14X 10-22 


6.77x10-75 


7.683 






14(3 


2i8Rn 


35.302 


2.41 X 10-23 


3.46x10-36 


26.374 






220 


210po 


32.841 


5.05X 10-23 


9.99x10-34 


41.279 








208pb 


21.528 


3.78X 10-22 


2.06x10-22 


62.309 






26Ne 


206pjg 


28.832 


2.98 X 10-24 


1.02x10-29 


57.966 
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Table 2: Continued 



Parent 
nucleus 


Emitted 
cluster 


Daughter 
nucleus 


Log Ti/2 
sec 


Preformation 
probability 

Po 


Decay 
constant 


Q value 
MeV 




234U 


^He 


230rpjj 


14.443 


1.64 xlO-"*^ 


2.50X 10-15 


4.859 






lOBc 




86.579 


1.04X 10-23 


1.83X 10-^^ 


6.713 










42.214 


9.38X 10-26 


4.24x10-^3 


24.514 






22Q 


212po 


38.903 


1.35X 10-25 


8.66x10-^° 


39.231 






24Ne 


210pb 


28.775 


7.96X 10-26 


1.16X 10-29 


58.825 






26Ne 


208 p]-) 


26.411 


1.14X 10-23 


2.69X 10-2^ 


59.465 




236 u 


^He 

lORe 

Uq 
22Q 

24Ne 

2«Ne 


232rpjj 

226Ra 
222Rn 

214pQ 
212pb 

2-'-"Pl) 


16.521 

93.871 
47.102 
42.807 
34.18 

31.117 


8.58X 10-1° 

1.06 X 10-23 
1.57x 10-26 
1.77X 10-26 
8.13X 10-28 

1.50 X 10-2'' 


2.09X 10-1^ 

9.33X 10-95 

6.74X 10-^^ 
1.08X 10-^3 
4.58X 10-35 

2.65 X 10-'''2 


4.574 
6.171 
23.054 
37.631 
55.944 
56.715 
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400 



CO 



CM 



o 




32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 

N (daughter) 



250 



Z=40 (daughter) 



(0 



o 



200- 



150- 



- 100- 



50- 







min 



I ' I ' I I I ' I I I ' I I I ' I ' I ' I ' I I I 
32 36 40 44 48 52 56 60 64 68 72 76 80 84 



N (daughter) 




N^(daughter) 



